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We propose a model for resonant Josephson tunneling through quantum dots that accounts for 
Cooper pair-breaking processes in the superconducting leads caused by a magnetic field or spin-flip 
scattering. The pair-breaking effect on the critical supercurrent Jc and the Josephson current- 
phase relation I{^p) is largely due to the modification of the spectrum of Andreev bound states 
below the reduced (Abrikosov-Gorkov) quasiparticle gap. For a quantum dot formed in a quasi-one- 
dimensional channel, both Jc and /(v?) can show a significant magnetic field dependence induced by 
pair breaking despite the suppression of the orbital magnetic field effect in the channel. This case 
is relevant to recent experiments on quantum dot Josephson junctions in carbon nanotubes. Pair- 
breaking processes are taken into account via the relation between the Andreev scattering matrix 
and the quasiclassical Green functions of the superconductors in the Usadel limit. 



PACS numbers: 74.50.-|-r,73.63.-b 



INTRODUCTION 



Since its discovery the Josephson efFeclii has been stud- 
ied for a variety of superconducting weak links^ii^iii. The 
research has recently entered a new phase with the ex- 
perimental realization of quantum dot weak links ex- 
ploiting electronic properties of finite-length carbon nan- 
otubes coupled to superconducting leads^i^ii. In particu- 
lar, for the first time since its theoretical prediction^i^ii^ 
resonant Josephson tunneling through discrete electronic 
states has been observed in carbon nanotube quantum 
dots^. As demonstrated in Refs. 00, the novel type of 
weak links exhibits transistor-like functionalities, e.g. a 
periodic modulation of the critical current with a gate 
voltage tuning successive energy levels in the dot on- 
and off- resonance with the Fermi energy in the leads. 
This property has already been implemented in a recently 
proposed carbon nanotube superconducting quantum in- 
terference device (CNT-SQUID)ii with possible applica- 
tions in the field of molecular magnetism. 

Motivated by the experiments on resonant Josephson 
tunneling, in this paper we investigate theoretically how 
robust it is with respect to pair-breaking perturbations in 
the superconducting leads. Cooper pair breaking can be 
induced by a number of factors, e.g. by paramagnetic im- 
purities^^, an external magnetic field^'^ or by structural 
inhomogeneities producing spatial fluctuations of the su- 
perconducting coupling constanti^. It can cause a drastic 
distortion of the Bardeen-Cooper-Schrieffer (BCS) super- 
conducting state, which manifests itself in the smearing 
of the BSC density of states leading to gapless supercon- 
ductivit}ii2.ii. 

While the pair-breaking effect on bulk superconductiv- 
ity is now well understood, its implications for quantum 
superconducting transport have been studied to a much 
lesser extent [see, e.g. Refs. 4,15,16] which to our knowl- 
edge does not cover Josephson tunneling through quan- 
tum dots. On the other hand, in low-dimensional systems 



pair-breaking effects may be observable in a common ex- 
perimental situation when, for instance, a carbon nan- 
otube weak link is subject to a magnetic field. Since the 
orbital field effect in the quasi-one-dimensional channel is 
strongly suppressed, pair breaking in the superconduct- 
ing leads can be the main source of the magnetic field 
dependence of the Josephson current. This situation is 
addressed in our work. 

The influence of pair breaking on the Josephson cur- 
rent can not, in general, be accounted for by mere sup- 
pression of the order parameter in the superconducting 
leads. As was pointed out in Ref. [H, it is a more sub- 
tle effect involving the modification of the spectrum of 
current carrying states in the junction, in particular, the 
subgap states usually referred to as Andreev bound states 
(ABS) ^°d'^ . We illustrate this idea for quantum dot junc- 
tions in the simple model of a short superconducting 
constriction with a scattering region containing a sin- 
gle Breit-Wigner resonance near the Fermi energy. The 
Josephson current is calculated using the normal-state 
scattering matrix of the system and the Andreev reflec- 
tion matri5*^ii^. Unlike Refs. [olfTol we focus on dirty su- 
perconductors for which the Andreev matrix can be quite 
generally expressed in terms of the quasiclassical Green 
functions^^, allowing us to treat pair breaking in the su- 
perconducting leads nonperturbatively. Although we ac- 
count for all energies (below and above the Abrikosov- 
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FIG. I: Scheme of a superconducting constriction with a nor- 
mal scattering region A*'. The arrows indicate the electrons 
(e) and holes (h) incident on and outgoing from A'^. 
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Gorkov gap Ag), it turns out that the behavior of the 
Josephson current can be well understood in terms of 
a pair-breaking- induced modification of the ABS, which 
depends sensitively on the relation between the Breit- 
Wigner resonance width T and the superconducting pair- 
ing energy A. Both the critical supercurrent and the 
Josephson current-phase relation are analyzed under ex- 
perimentally realizable conditions. 



II. MODEL AND FORMALISM 

We consider a junction between two superconduc- 
tors Si and ^2 adiabatically narrowing into quasi-one- 
dimensional ballistic wires S'l and 5*2 coupled to a nor- 
mal conductor N [Fig. [T]. The transformation from the 
superconducting electron spectrum to the normal-metal 
one is assumed to take place at the boundaries SiS'i 
and S2S2, implying the pairing potential of the form-: 
A{x) = Ae'^i for x < ~L/2, A{x) = for \x\ < L/2 
and A(x) = Ae*'^^ for x > L/2 with the order parameter 
phase difference ip = (p2 — ipi and the junction length 
L <C hvp/A [vp is the Fermi velocity in Si^2)- 

The Josephson coupling can be interpreted in terms 
of the Andreev process^ whereby an electron is retro- 
reflected as a Fermi-sea hole from one of the super- 
conductors with the subsequent hole-to-electron conver- 
sion in the other one. Such an Andreev reflection cir- 
cle facilitates a Cooper pair transfer between S'l and 5*2. 
Normal backscattering from disordered superconducting 
bulk into a single-channel junction is suppressed due to 
the smallness of the junction width compared to the elas- 
tic mean free path i. The N region in the middle of the 
junction is thus supposed to be the only source of nor- 
mal scattering. In such type of weak links the Josephson 
current is conveniently described by the scattering ma- 
trix expression of Refs.[Tolf20l that can be written at finite 
temperature T as the following sum over the Matsubara 
frequencies aj„ = {2n + \)'KkBT [Ref . [20| : 
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Here s n (E) is a 4 x 4 unitary matrix relating the incident 
electron and hole waves on the N region to the outgoing 
ones [Fig. [1]. It is diagonal in the electron-hole space: 



sn 



See{E) 
Shh{E) 



AE) 



rii{E) ti2{E) 
t2i{E) r22{E) 



The matrix See {E) describes electron scattering in terms 
of the reflection and transmission amplitudes, r.jk{E) 
and tjk{E), for a transition from 5^ to S'j {j,k — 1,2). 
The hole scattering matrix is related to the electron one 
by Shh{E) = s*g(— -B). The Andreev scattering matrix 
sa{E) is ofl^-diagonal in the electron-hole space: 



SA 





She{E) 



eh{E) 





(2) 



where the 2x2 matrices She{E) and Seh{E) govern the 
electron-to-hole and hole-to-electron scattering off the su- 
perconductors. Equation ([T|) is valid for all energies as 
long as normal scattering from the superconductors is 
absenti^i^. 

In Ref. [l^ the Andreev matrix ^ was obtained by 
matching the solutions of the Bogolubov-de Gennes equa- 
tions in the wires S[ 2 to the corresponding solutions in 
impurity-free leads. Gorkov's Green function formalism 
in combination with the quasiclassical theory^ allows 
one to generalize the results of Ref. [l3 to dirty leads with 
a short mean free path £ ^ Hvp/A. In the latter case the 
matrices She{E) and Seh{E) can be expressed in terms of 
the quasiclassical Green functions of the superconductors 
as foUowsi^: 



Seh 



ME) 

gi(E)+i 






/2(g) 

S2(B) + ] 



5 She 



-fliE) 

gi{E)+i 




~fUE) 
92(£) + l 



Here gi^2 and /i,2 (/i 2) ^^^i respectively, the normal and 
anomalous retarded Green functions in 51^2- These ma- 
trices are diagonal in the electrode space due to a local 
character of Andreev reflection in our geometry. 

Neglecting the influence of the narrow weak link on the 
bulk superconductivity, we can use the Green functions 
of the uncoupled superconductors 51,2 described by the 
position-independent Usadel equation^, 



ih 

Tpb 



Efs + Aj + — T3gjT3 , gj 



0, 



(3) 



with the normalization condition g| = fg for the matrix 
Green function 



9j 



9] fj 



Ac'"^:'- 
-Ae-''^^ 



1,2. 



Here tq and are the unity and Pauli matrices, re- 
spectively, and [...,...] denotes a commutator. Equa- 
tion ([3|) accounts for a finite pair-breaking rate 
whose microscopic expression depends on the nature of 
the pair-breaking mechanism. For instance, for thin su- 
perconducting films in a parallel magnetic field, = 

(i;F^/18)(7rdB/$o)^ [Ref.lll where d is the film thickness 
and is the flux quantum. For paramagnetic impuri- 
ties, Tph coincides with the spin-flip time^^i. In the case of 
the spatial fluctuations of the superconducting coupling, 
is proportional to the variance of the fluctuationsii. 
From Eq. ([3]) one obtains the Green functions 



9} = 



= "e-'^^/„ /j--e-2»'^V„ (4) 



-r ^ u\\ , 

A V yr^i? 



c 



(5) 



where, following Refs. [T^fTsl . we introduce a dimension- 
less pair-breaking parameter C,. The matrices and s-^f^ 
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can be expressed using Eqs. Q as follows: 

















Seh = a 





eiv2 










= u — \/v? — 1. 



(6) 
(7) 



We note that pair breaking modifies the energy depen- 
dence of the Andreev reflection amplitude a according to 
the non-BCS Green functions ([4]) and ([5]). A few words 
concerning the applicability of this result are due here. 

First of all, there is no restriction on energy E, e.g. for 
C < 1, equations ^ and ([7]) are valid both below and 
above the reduced (Abrikosov-Gorkov) quasiparticle gap 

Ag = A (l - ^2/3)^/2_ In particular, for \E\ < Ag one 
can show that u is real and |u| < (1 — < l 

[Ref. ISj], corresponding to perfect Andreev reflection 
with a = exp(— i arccos(u)). Since in the Usadel limit 
£ <C vpTpb, normal scattering from the superconduc- 
tors is suppressed due to the smallness of the junction 
width also in the presence of pair breaking. The ab- 
sence of normal transmission at \E\ < Ag is consistent 
with the Abrikosov-Gorkov approach assuming no impu- 
rity states inside the gap and the validity of the Born 



approximatio] 
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For \E\ > As the relevant solution 



of Eq. ^ is complex and has positive Imu related to the 
density of states of the superconductor^"^. Equations ([5]) 
and ([7]) are thus the generalization of the known result 
a = (E/Aq) - ■s/{E/AoY - 1 [Ref. 22] for transparent 
point contact, where Ao = A|^=o is the ECS gap. It 
is convenient to measure all energies in units of Aq for 
which equations (U)-© should be complemented with 
the self-consistency equation for A. At T = 0, the case 
we are eventually interested in, this equation can be writ- 
ten as— li^: 



ln(Co/C) 
ln(Co/C) 



-7rC/4, 



C< 1, 

1/(20 -ln(C - 



v/C^/(2C)-ln(C+VC^ 
(C/2)arctan(l/VC^) , 



(8) 

1)- (9) 
C> 1, 



with C being now a function of a new pair-breaking pa- 
rameter Co — ^/(''pbAo) ranging from zero to the critical 
value Co = 0.5 at which C, — co and A — Oi^iH. 

Inserting Eqs. ^ and ^ for sa{E) into Eq. ^ and 
taking the limit T — + we obtain the Joscphson current 
for an arbitrary sn{E) as 



I = -- 



4e 



d 

du! In { 1 

dtp 



Det SeeiE) Det s*ee{-E) - 
[niiE)rl^i-E) + r22iE)r;^{-E) + 
e-'n2i{E)tU-E) + e^n,2{E)tl^{-E)\ ]^ 



(10) 
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FIG. 2: Phase dependence of the Andreev bound state for a 
broad resonant level with V = 15Ao close to the Fermi energy 
(-Br = O.ir); dashed line shows the normalized gap for a given 
value of the pair-breaking parameter Co. 



III. ANDREEV BOUND STATES IN A 
RESONANT JUNCTION 

Let us assume that the N region is a small quantum 
dot and electrons can only tunnel via one of its levels 
characterized by its position Er with respect to the Fermi 
level and broadening F. For the simplest Breit-Wigner 
scattering matrix with rii ~ r22 = {E—Er)/{E—Er+iT) 
and ti2 ~ t2i ~ r/i(i? — E,. + iT), equation pU]) reads 



/ = — (2e//i)T sill (fi J dujx 





(11) 



7^ + r 1 



r/A 



-l + Tsin^f^) 



E—iuj. 



where T = 1 - 7^ = rVl^r + T^) is the Breit-Wigner 
transmission probability at the Fermi level. The param- 
eter F/A accounts for the energy dependence of the res- 
onant superconducting tunneling. In Eq. (jllD the inte- 
grand has, in general, poles given by the equation 



7^ + r 1 



r/A 



l-rsin'(f).(12) 



Along with Eq. ((S]) they determine the energies of the An- 
dreev bound states (ABS) localized below the Abrikosov- 
Gorkov gap Ag = A (l — C^/'^)'^^^. It is instructive to un- 
derstand how the pair breaking modifies the ABS spec- 
trum since this is reflected on both the current-phase 
relation I{tp) and the critical current Ic = max /(</?). 

We start our analysis with an analytically accessible 
case of an infinitely broad resonant level, F/A — > oo, 
where Eq. reduces to = 1 — T siv?{ip/2), yielding 
the ABS energies ±E{ip) [see, Eq. Q]: 



E{lp) = a Jl - Tsin\Lp/2) 



1 



c 



T|sin(v,/2)| 



(13) 



Requiring £'(1^9) < Ag we find that the ABS exist in 
the phase interval where sin^ (1^9/2) > C,'^/'^ /T and only 
if C^/^ < T. The numerical solution of Eqs. ©, ® 
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FIG. 3: Phase dependence of the Andreev bound state for a 
narrow resonant level with F — 0.3Ao and Er — O.IF. 
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and confirms that the interval of the existence of 
ABS gradually shrinks from < (/J < 27r to a narrower 
one with increasing pair breaking [see, Fig. [2] . Outside 
this interval the Josephson current is carried by the con- 
tinuum states {E > Ag) alone, which is automatically 
accounted for by Eq. (fTTj) . An equation of the same form 
as Eq. (jl3p was derived earlier for a nonresonant system 
and by a different method-'^^. 

By contrast, the ABS spectrum for a narrow reso- 
nant level turns out to be much less sensitive to pair 
breaking. Indeed, under condition F/A ^ 1 — C equa- 
tions ^ and (fT2|) reproduce the known result, E{(p) — 

V^£;2 -f r2y^l -rsin2((^/2) [Refs. ii]. In particular, 
for Er ^ the ABS exist within the resonance width, 
E[ip) < F and are separated from the continuum by a 
gap Ag - F. Solving Eqs. (O, ^ and p2|l numerically, 
we find that until this gap closes at a certain value of 
the ABS spectrum remains virtually intact [see. Fig. [3]. 
For bigger (q, the spectrum gets modified in a way simi- 
lar to the previous case [cf., third panels in Figs. [2] and [3]. 
In the case of a very narrow resonance, the characteristic 
value of Co is « 0.45, corresponding to C ~ 1, i-e. to the 
onset of gapless superconductivit y^^i^^ . 



FIG. 4: On-resonance critical current vs. pair-breaking pa- 
rameter (^0 = h/{TpbAo) for different F/Ao. The behavior of 
the normalized order paramete r^^'^^ is shown, for comparison, 
in red. 



Using Eqs. (O, ([9]) and we are able to analyze the 
critical current Ic = max I{ip) in the whole range of the 
pair-breaking parameter, < Co ^ 0-5 [see. Fig. 0]. In 
line with the discussed behavior of the Andreev bound 
states, for a narrow resonance, F/Aq <C 1, the critical 
current starts to drop significantly only upon entering 
the gapless superconductivity regime 0.45 < Co < 0.5. 
On the other hand, for a broad resonance, T/Aq ^ 
1, the suppression of Ic is almost linear in the whole 
range. We note that in both cases the behavior of Ic 
strongly deviates from that of the bulk order parameter 
(red curve)J^-'" - largely due to the pair-breaking effect 
on the ABS. In practice, the /c(Co) dependence can be 
measured by applying a magnetic field [the case where 
Co = {B/B^f and = (^o/nd).J\^K^Jhv^i] in an ex- 
periment similar to Ref. [y where a quantum dot, defined 
in a single-wall carbon nanotube, was strongly coupled 
to the leads with the ratio F/Aq ~ 10. Carbon nanotube 



IV. CRITICAL CURRENT AND 
CURRENT-PHASE RELATION: RESULTS AND 
DISCUSSION 

For numerical evaluation of the Josephson current (fTT|) 
we first put E = iLO in Eq. ((5]) and then make the transfor- 
mation u — > iv, yielding cj/A — — C/vT+i^). Using 
this relation, in Eq. (fTTj) we change to the integration 



over V with the Jakobian dhj/dv = A[l — C/(l + i'" , 

CO , V 

/=(2eA//i)rsin(p/dz. (l- (^^^j X (14) 



l-Tsin^d) 



Positiveness of uj in Eq. p4|) enforces the choice of the 
lower integration limit: i^o = for C ^ 1 ^-nd vt^ = 
.JC - 1 for C > 1- 
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FIG. 5: On-resonance current-phase relation for different 
values of the pair-breaking parameter Co and F = Aq. 
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FIG. 6: Critical current vs. resonant level position: ~ 
(black), Co = 0.25 (blue) and Co = 0.45 (red). 

quantum dots with lower T/Aq values are accessible ex- 
perimentally, too^iii. 

We also found that the crossover between the gapped 
and gapless regimes is accompanied by a qualitative 
change in the shape of the Josephson current-phase rela- 
tion /((/?) as demonstrated in Fig. [5] for the on-resonance 
case Er = and F = Aq. The I{(p) relation is anhar- 
monic as long as the junction with Ag ^ supports 
the ABS (black and blue curves). The vanishing of the 
ABS upon entering the gapless regime leads to a nearly 
sinusoidal current-phase relation (red curve). A closely 
related effect is demonstrated in Fig. [5] showing the mod- 
ification of the critical current resonance lineshape with 
the increasing pair-breaking strength. In the absence of 
pair breaking it is nonanalytic near Er = (black curves) 
reflecting the anharmonic I{ip) due to the ABS in a trans- 



parent channel^'W. On approaching the gapless regime 
this singularity is smeared out (red curves), which is ac- 
companied by the suppression of the Ic amplitude. At fi- 
nite temperatures T <^ A/ks the pair-breaking-induced 
smearing of the resonance peak will enhance the usual 
temperature effect. 

In conclusion, we have proposed a model describing 
resonant Josephson tunneling through a quantum dot 
beyond the conventional BCS picture of the supercon- 
ducting state in the leads. It allows for nonperturbative 
treatment of pair-breaking processes induced by a mag- 
netic field or paramagnetic impurities in diffusive super- 
conductors. We considered no Coulomb blockade effects, 
assuming small charging energy in the dot Ec <C Aq, F, 
which was, for instance, the case in the experiment of 
Ref . [el. Our predictions, however, should be qualitatively 
correct also for weakly coupled dots with F < Eq ^ Aq 
at least as far as the dependence of the ctitical super- 
current on the pair-breaking parameter is concerned. In- 
deed, for a narrow resonance the Andreev bound states 
begin to respond to pair breaking only when the gap Ag 
becomes sufficiently small [see. Fig. |3] so that for a fi- 
nite Ec ^ Aq one can expect a sharp transition to the 
resistive state, too, similar to that shown in Fig. [1] for 
r/Ao « 1. 
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